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In  this  paper  the  linear  and  nonlinear  theory  of  the  electron  cyclotron  maser  instability  is  con- 
sidered. The  configuration  used  to  study  the  maser  instability  consists  of  relativistic  electrons 
gyrating  about  and  drifting  along  a uniform  magnetic  field  within  a parallel  plate  waveguide.  Rela- 
tivistic effects  associated  with  the  gyrating  electrons  are  responsible  for  excitation  of  the  transverse 
electric  mode  in  the  waveguide.  Linear  theory  shows  that  the  growth  rate  maximizes  when  the  actual 
beam  velocity  coincides  with  the  axial  wave  group  velocity  of  the  excited  electromagnetic  wave.  This 
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allows  us  to  perform  thelionlinear  analysis  in  a frame  where  both  the  axial  wave  number  and  axial 
beam  velocity  vanish.  ^We  have  found  that  the  maser  instability  exists  only  if  the  perpendicular  beam 
energy  exceeds  a threshold" value.  Our  analysis  also  describes  the  temporal  nonlinear  evolution  of  the 
field  amplitude  and  frequency  of  a single  excited  wav&.  The  nonlinear  wave  dynamics  are  self  con- 
sistently determined  from  the  nonlinear  particle  orbits  through  the  force  and  wave  equations.'JThe 
nonlinear  analysis  shows  that  there  are  two  possible  mechanisms  for  the  saturation  of  the  unstable 
wave: 

i.  depletion  of  the  available  free  energy  associated  with  the  rotating  particles 

ii.  phase  trapping  of  the  gyrating  electrons  in  the  wave. 

The  initial  beam  parameters  determine  which  of  the  two  mechanisms  is  responsible  for  saturation. 
Competition  between  the  two  saturation  mechanisms  leads  to  a peaking  in  the  energy  conversion 
efficiency  as  a function  of  beam  energy.  Numerical  results  of  the  nonlinear  formalism  show  that 
energy  conversion  efficiencies  from  the  particles  to  the  wave  can  be  as  high  as  60  percent  in  the  beam 
frame. furthermore,  by  appropriately  contouring  the  external  magnetic  field,  among  other  things, 
efficiencies  as  high  as  70  percent  can  be  realized. 
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THE  LINEAR  AND  SELF  CONSISTENT  NONLINEAR  THEORY 
OF  THE  ELECTRON  CYCLOTRON  MASER  INSTABILITY 

Section  I.  Introduction 

In  the  cyclotron  maser  instability  electromagnetic  radiation  is 
generated  by  relativistic  electrons  gyrating  about  an  external  magnetic 
field.  The  bunching  process  in  the  maser  instability  is  a relativistic 
effect  due  to  the  energy  dependent  electron  cyclotron  frequency.  The 
excited  wave  frequency  is  near  the  Doppler  shifted  electron  cyclotron 
frequency  of  the  gyrating  electrons.  The  free  energy  for  the  maser  in- 
stability resides  in  the  rotational  energy  of  the  electrons.  To  make 
this  source  of  energy  from  a streaming  electron  beam  available  for  the 
generation  of  radiation,  a redistribution  of  electron  energy  is  needed. 
For  example,  converting  streaming  to  rotational  energy  can  be  accom- 
plished by  propagating  the  electron  beam  through  a spatially  varying 
external  magnetic  field.  This  conversion  process  is  the  classical 
equivalent  of  producing  inverted  population  levels  in  molecular  masers 
by  pumping. 

The  electron  cyclotron  maser  instability  is  particularly  promising 
for  efficient  generation  of  millimeter  and  submillimeter  radiation.  In 
this  part  of  the  electromagnetic  spectrum  more  conventional  devices  are 
severely  limited  in  both  power  and  efficiency.  Such  devices  as  the 
traveling  wave  tube  or  magnetron  rely  on  a slow  wave  structure  for  their 
operation.  The  radiation  wavelength  from  these  devices  is  typically  on 
the  order  of  the  structure  period.  Power  density  as  well  as  voltage 
breakdown  considerations  place  a lower  limit  on  the  dimensions  of  the 
structure  depending  on  the  efficiency  and  radiated  power  requirements. 
Therefore,  these  slow  wave  devices  are  not  particularly  suited  for 
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efficient  generation  of  large  fluxes  of  radiation  at  millimeter  and  sub- 
millimeter wavelengths.  The  electron  cyclotron  maser  mechanism,  on  the 
other  hand,  does  not  rely  on  the  fine  structure  of  a waveguide  or  cavity 
and  thus  efficient  operation  at  millimeter  and  submillimeter  wavelengths 
is  possible. 

The  linear  mechanism  of  the  cyclotron  maser  was  first  proposed  by 
Twiss1  classically  and  later  quantum  mechanically  by  Schneider.2  Shortly 
after  the  work  of  Twiss,  Gapanov  published  a paper  on  the  classical 
theory  of  the  cyclotron  maser.9  In  1959,  Pantell  published  what  was 
perhaps  the  first  experimental  work  involving  the  electron  cyclotron 
maser  mechanism.4  In  this  experiment  radiation  at  2.5  to  4.0  Ghz  was 
produced  from  a 1 kv,  3 pP  electron  beam.  Pantell  suggested  that  the 
radiation  was  caused  by  the  coupling  of  the  TE^  waveguide  mode  to  the 
backward  traveling  cyclotron  wave  on  the  magnetized  beam.  It  was  Hirsh- 
field  and  Wachtel  who  performed  the  first  experiment  that  definitely 
confirmed  the  existence  of  the  electron  cyclotron  maser  mechanism. 5 
Apparently,  because  of  a lack  of  understanding  of  the  nonlinear  dynamics, 
the  observed  efficiencies  in  the  initial  experiments  were  low,  percent. 
In  this  early  work  the  cyclotron  maser  was  operated  both  as  an  oscillator 
and  as  an  amplifier  at  millimeter  wavelengths. 5-9  The  radiated  power 
achieved  in  these  experiments  was  as  high  as  one  kilowatt.8  Nonpulsed 
low  current  electron  beams  were  used  throughout  the  early  experiments. 

The  advent  of  intense  pulsed  relativistic  electron  beams  has  re- 
newed interest  in  the  cyclotron  maser  instability  as  a source  of  high 
power  microwave  radiation.  These  intense  electron  beams  use  a field 
emission  diode  and  typically  carry  a current  of  100  kA,  a voltage  of 
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■sf 5 MeV  and  last  for  approximately  50  nsec.  Microwave  radiation  at 
multi-megawatt  levels10  and  as  high  as  one  gigawatt11  has  been  produced 
using  these  intense  relativistic  electron  beams.  By  operating  the  cyclo- 
tron maser  as  a single-stage  amplifier,  a 17  dB  gain  was  obtained.12 
In  general  the  emitted  radiation  was  narrow  band  with  a coherence  time 
> 50  nsec,  and  was  also  spatially  coherent.13  None  of  the  experiments 
using  pulsed  intense  electron  beams  attained  efficiencies  greater  than 
about  2 percent.  Our  analysis  shows  that  there  are  a number  of  reasons 
for  the  observed  low  efficiencies  with  intense  relativistic  electron 
beams . 

The  highest  reported  efficiencies  with  a cyclotron  maser  device 
were  obtained  in  the  Soviet  Union  using  thermionic  diodes.14*15  In  the 
experiments  by  Zaytsev  et  al. ,14  12  kw  cw  at  2.78  mm  was  reported  using 
the  fundamental  cyclotron  mode.  Experimental  efficiencies  in  this  case 
were  31  percent.  Interaction  at  the  second  cyclotron  harmonic  produced 
an  output  power  of  1.5  kw  cw  at  .92  mm  with  an  efficiency  of  6.2  percent. 
Using  the  second  cyclotron  harmonic  Kisel  et  al.15  obtained  output  powers 
of  10  kw  cw  at  8.9  mm  with  efficiencies  of  40  percent.  In  the  pulsed 
electron  mode,  power  levels  of  30  kw  at  43  percent  efficiency  were 
reported. 

Work  on  the  linear  theory  of  the  cyclotron  maser  shows  that  the 
instability  is  due  to  the  coupling  of  the  supraluminous  TE  waveguide 
mode  and  a beam  cyclotron  mode.10*15-13  The  instability  may  be  either 
absolute  or  convective  in  nature.  For  a review  of  the  theory  on  excited 
nonlinear  oscillators  as  applied  to  the  cyclotron  maser,  see  the  article 
by  Gapanov  et  al.20 
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To  our  knowledge,  Che  first  nonlinear  evaluation  of  the  electron 
cyclotron  maser  mechanism  involved  numerical  integration  of  the  elec- 
tron orbits  in  fields  of  either  constant  amplitude  and/or  constant 
f requency . 1 f > 2 5 These  nonlinear  theories  do  not  fully  treat  the 
particles  and  wave  dynamics  in  an  inherently  self-consistent  manner. 

In  this  paper  a comprehensive  study  of  the  self-consistent  non- 
linear evolution  of  the  electron  cyclotron  maser  instability  is  pre- 
sented. We  have  included  both  a time-dependent  frequency  shift  and  a 
time  dependent  field  amplitude  in  our  analysis.  A rather  simple  physical 
picture  of  the  nonlinear  behavior  of  this  instability  shows  that  there 
are  two  mechanisms  responsible  for  wave  saturation.  The  analysis  is 
self  consistent  in  that  the  particle  and  wave  dynamics  are  treated  as  a 
unit.  Our  parameter  study  indicates  that  the  maser  instability  can  be 
an  efficient  mechanism  for  the  production  of  high  power  radiation, 
particularly  at  millimeter  and  submillimeter  wavelengths.  Saturation 
efficiencies  at  the  fundamental  waveguide  mode  and  cyclotron  mode  are 
obtained.  A method  of  externally  controlling  the  saturation  process, 
and  thereby  further  improving  conversion  efficiencies,  is  discussed. 

The  physical  mechanism  of  the  electron  cyclotron  maser  is  described 
in  Section  II.  Properties  of  the  linear  dispersion  relation  are  dis- 
cussed in  Section  III.  Here  simple  expressions  for  the  linear  growth 
rate  and  frequency  shift  are  obtained.  The  linear  theory  indicates 
that  the  growth  rate  maximizes  when  the  axial  beam  velocity  is  equal  to 
the  group  velocity  of  the  excited  wave.  This  fact  can  be  used  to 
simplify  the  analysis  by  performing  all  calculations  in  the  beam  frame 
which  now  coincides  with  the  cutoff  frame  where  the  axial  wave  number, 
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k , vanishes.  The  results  are  then  transformed  back  to  the  laboratory 
o 

frame.  In  Section  IV  we  derive  the  equations  for  the  nonlinear  growth 
rate  and  frequency  shift  of  the  excited  wave.  The  nonlinear  single  wave 
formalism  treats  both  the  particle  and  wave  dynamics  self  consistently. 

The  equivalence  of  our  temporal  model  with  the  steady  state  spatial  growth 
of  a wave  in  an  amplifier  is  pointed  out.  A number  of  simplifying  assump- 
tions, such  as  the  single  wave  approximation  and  neglect  of  initial  beam 
temperature  and  space  charge  effects,  have  been  made  in  our  analysis.  The 
conditions  for  the  validity  of  these  approximations  are  discussed  in 
Section  V.  Section  VI  describes  the  two  nonlinear  saturation  processes 
of  the  cyclotron  maser  instability.  There  we  show  that  the  depletion  of 
free  energy  and/or  particle  phase  trapping  is  responsible  for  wave  satura- 
tion. The  dominant  mechanism  depends  on  the  initial  choice  of  beam  param- 
eters. Simple  analytic  expressions  for  the  field  amplitude  at  saturation 
are  obtained.  From  the  slow  time  scale  orbit  equations  we  obtain  a con- 
stant of  the  motion  for  a fixed  amplitude  and  frequency  field.  This 
constant  is  used  to  examine  the  nonlinear  particle  orbits  in  phase  space. 
Section  VII  contains  the  numerical  results  of  our  formalism  for  a wide 
range  of  beam  parameters.  These  results  include  field  amplitudes  and 
efficiencies  at  saturation.  Comparison  of  these  results  with  the  analy- 
tic expressions  given  in  Section  VI  is  made.  Methods  for  improving  the 
efficiency  at  saturation  are  also  discussed  and  examples  given.  Section 
VIII  contains  a brief  discussion  of  our  results  as  well  as  possible 
applications  for  the  cyclotron  maser.  The  low  experimental  efficiencies 
with  intense  relativistic  electron  beams  for  the  cyclotron  maser  are 
discussed  in  the  light  of  our  theoretical  results.  A detailed  example 
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showing  how  our  results  can  be  applied  to  practical  situations  to  esti- 
mate power  efficiencies  and  levels  of  the  radiation  field  is  given.  In 
Appendix  A we  show  that  the  linear  dispersion  relation  can  be  recovered 
from  the  nonlinear  formalism.  Because  our  analysis  is  performed  in  the 
beam-cutoff  frame,  we  include  a derivation  of  the  efficiency  transforma- 
tion from  the  beam  frame  to  the  laboratory  frame  in  Appendix  B. 
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Section  II.  Physical  Mechanism 


Insight  into  the  physical  mechanism  responsible  for  the  electron 
cyclotron  maser  instability  can  be  obtained  from  the  particle  trajec- 
tories shown  in  Figs.  1(a)  and  1(b).  These  figures  show  the  orbit 
dynamics  of  sample  electrons  initially  uniformly  distributed  along  a 
gyro  orbit.  The  electrons  are  rotating  in  the  counter  clockwise  direc- 
tion about  a uniform  and  constant  magnetic  field,  B e . Without  loss 

o z 

of  generality  we  assume  the  electron  velocity  in  the  z direction  to  be 

zero.  The  initial  radius  of  the  sample  electron  ring  is  the  Larmor  radius 

r T = v / (Q  /v  ) where  v is  the  initial  perpendicular  velocity, 
oL  ox  o ox  ox  J 

= ie!Bo/moC  is  the  nonrelativistic  cyclotron  frequency  and 

Y = (1  - v2  /c2)"1/2  is  the  relativistic  factor, 
ox  oi 

The  introduction  of  a small  constant  amplitude  electric  field, 

Ey(t)  = cos  (wot)e  , as  shown  in  Fig.  1(a),  will  alter  the  particle 
orbits.  We  first  examine  the  particle  trajectories  when  the  frequency 
of  the  electric  field  is  equal  to  the  initial  relativistic  cyclotron 
frequency,  = fi  /y  . The  time  rate  of  change  of  the  particle  energy 
is  de^/dt  = - ]e|  v^(t)E  (t)  where  v^(t)  is  the  y component  of  particle 
velocity.  With  the  initial  choice  of  field  direction  shown  in  Fig.  1(a), 
particles  8,  1 and  2 will  lose  energy  and  tend  to  spiral  inward.  The 
relativistic  cyclotron  frequency  of  these  particles  will  increase,  since 
Y,  decreases,  and  their  phase  will  tend  to  slip  ahead  of  the  wave. 
Particles  4,  5 and  6,  on  the  other  hand,  will  gain  energy,  their  cyclo- 
tron frequency  will  decrease  and  they  will  tend  to  spiral  outward.  The 
phase  of  these  particles  will  tend  to  slip  behind  the  wave.  After  an 
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integral  number  of  wave  periods,  the  particles  will  become  bunched  around 
the  positive  y-axis.  Particle  bunching  is,  therefore,  caused  by  rela- 
tivistic effects,  since  the  rotational  frequency  of  the  electrons  is 
energy  dependent. 

To  obtain  a net  exchange  of  energy  between  the  particles  and  wave, 

u must  be  slightly  greater  than  ft  /y  . If  u £;ft  /y  , the  particles 
o o 01  o o OJ. 

on  the  average  traverse  a coordinate  space  angle  less  than  2it  in  a wave 

period,  2 tt/oj  . All  the  particles  will  then  slip  behind  the  wave  and  the 
o 

distribution  of  particles  after  an  integral  number  of  wave  periods  will 

be  bunched  in  the  upper  half  plane  as  shown  in  Fig.  1(b).  As  a result 

of  the  phase  slippage  between  the  particles  and  the  field,  the  net 

kinetic  energy  of  the  ensemble  of  particles  decreases.  Conservation  of 

total  energy  implies  that  the  field  amplitude  increases  resulting  in 

instability.  If  oo  remains  greater  than  ft  /v  the  particles  will  con- 
o o J. 

linue  to  slip  behind  the  wave. 

Depending  on  the  initial  beam  parameters,  the  group  of  bunched  beam 
particles  may:  i)  continue  to  slip  behind  the  wave,  or  ii)  initially 

slip  behind  the  wave,  reverse  itself  and  begin  to  oscillate  about  the 
positive  y-axis.  In  either  case,  the  bunched  particles  will  eventually 
appear  in  the  lower  half  of  the  transverse  plane  after  an  integral  number 
of  wave  periods.  When  this  occurs,  the  particles  will  gain  energy  and 
the  wave  amplitude  will  begin  to  decrease. 
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Section  III.  Linear  Theory 


Using  the  Maxwell-Vlasov  equations,  the  linear  theory  of  the  elec- 
tron cyclotron  uaser  instability  has  been  derived  in  both  cylindrical10 
and  cartesian  geometry. ^ The  electron  beam  and  waveguide  configuration 
employed  in  Ref.  19  and  shown  in  Fig.  2 is  basically  the  same  model  used 
in  the  present  nonlinear  analysis.  In  Fig.  2 a beam  of  gyrating  and 
drifting  particles  is  assumed  to  have  the  same  perpendicular  velocity, 
v^,  and  parallel  velocity,  vq(|  , with  respect  to  the  applied  axial 
magnetic  field  B^.  The  guiding  centers  of  the  particles  lie  midway  be- 
tween the  two  conducting  plates.  The  system  of  particles  and  fields  is 
assumed  independent  of  the  spatial  y variable.  The  field  components 

within  the  waveguide  are  those  of  a TE  mode  (i.e.,  E , B and  B are 

on  y z x 

the  only  nonvanishing  field  components).  The  functional  dependence  of 

the  fields  on  z and  t has  the  form  exp  (i(k  z - ait)).  The  linear  dis- 

o 

persion  relation  in  the  laboratory  frame  for  the  above  configuration  is 

w2  - c*  (k2  + k2)  = 6 „(w(-/y  ). 

o n nS,  b o 


(w 


k v ) Q . 
o o ii  n l 


B2  (u)2  - 


Oi 


c2k2) 

o 


(u)  - k v 

O O II 


- lW 


(i 


k v 

O O I! 


- ZQ  h )z 
o o 


(1) 


where  k = rn/2a  is  the  perpendicular  wave  number,  n = 1,  2,  3,  . . . 
n 

is  the  waveguide  mode  number,  5 . = (1  + (-  l)n+^) , w2  * An  lei  a /m  a 

n*.  b o o 

is  the  modified  plasma  frequency,  oq  is  the  surface  charge  density  of 

the  beam  in  the  x-y  plane,  i = 1,  2,  3 is  the  magnetic  harmonic  number, 

Y * (1  - v^  /c2  - v2  /c2)~!/“,  n = | e | B /m  c,  B = v /c, 
o on  oi  o ' 1 o o ox  ox 

Q , - x (l2/x2  - 1)  d(J. (x  ) )2/dx,  W - (dJ. (x  ) / dx  ) 2 , 
nx  n n x,  n nx  »■  n n 
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x =6  ck  / (ft  /y  ) and  J is  the  Bessel  function  of  order  H.  If  n + i 
n oj.  n o o l 

is  odd,  the  TE  mode  is  stable  since  the  right  hand  side  of  Eq . (1) 
on 

vanishes.--6  The  first  term  in  the  bracket  on  the  right  hand  side  of  (1) 

is  always  stabilizing  while  the  second  term  is  always  destabilizing. 

Furthermore,  the  frequency  of  the  unstable  wave  in  the  rest  frame  of  the 

electrons  is  slightly  greater  than  the  cyclotron  frequency, 

v (uj  - k v „ ) i ft  . We  have  seen  in  Section  II  that  this  requirement 
'o  o on  o 

is  necessary  for  a phase  slippage  to  occur  between  the  particles  and 

field.  It  is  clear  from  the  dispersion  relationship  in  Eq.  (1)  that  the 

maximum  growth  rate  of  the  modified  TE  mode  occurs  for  frequencies  and 

wave  numbers  near  the  intersection  of  the  vacuum  waveguide  mode 

w2  - c2(k2  + k2)  = 0,  and  the  cyclotron  mode,  uj  - k v - 1ft  /y  =0. 
on  o on  o o 

When  the  perpendicular  velocity  of  the  particles,  VQ1  » vanishes  the 
cyclotron  wave  is  a positive  energy  wave.  Since  the  waveguide  mode  is 
also  a positive  energy  wave,  there  can  be  no  instability  for  v =0. 

o J. 

If,  however,  v £ q,  the  cyclotron  mode  splits  into  a positive  and  nega- 
tive energy  polarity  wave  as  can  be  seen  from  the  fact  that 

(w  - k v - IS2  /y  )2  appears  in  the  denominator  of  the  destabilizing 
oouoo 

term  of  Eq.  (1).  The  coupling  of  the  negative  energy  cyclotron  wave 

and  positive  energy  waveguide  mode  is  responsible  for  the  instability. 

In  the  limit  of  vanishing  beam  density,  u -*•  0,  or  vanishing  perpen- 

b 

dicular  velocity,  v^  -*•  0,  the  opposite  energy  polarity  cyclotron  waves 
coalesce  and  become  degenerate. 

Seeking  a solution  to  Eq.  (1)  near  the  intersection  of  the  modes 

we  set  a a u + 6ui  where  1 I <<  u and  oj  = c(k2  + k2)1^  * 
o o ' o o o on 
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k v + HI  /y  . Substituting  oi  = o>  + <5u>  into  (1)  and  keeping  terms  to 

OO'IOO  oo 

order  (6uj  ) 5 yields  the  following  relation 


- 3 A2ia  6ui  + 3 Z A^ui  = 0 (2) 

o o o nl  o 

where 


^ u)2/(6y  w ))Q  ) 1/2  and  Z = (x  /i)2  (W  /Q  ) UR  /y  ) . 

o o o nx.  b o o nx  n£  n nx.ni  o o 

The  solutions  of  Eq.  (2)  yield  complex  roots  when  Z „ > (2/3)Au>  . 

nJc.  o 


This  condition  can  be  stated  as  a threshold  condition  for  instability, 


requiring  that  the  perpendicular  velocity  be  greater  than  a critical 
value , 2 7 


B >6 

ojl  i,crit 


where 


1/2  r 


11/4 


6 . - UR  /Y  )(Q  „/ W .)  1(2/27)6  .w2Q  ./(y  u injy  )|  /ck 

j.,crit  o o nil  nS.  L n£  b nil  o o o o J n 

At  the  critical  value  of  8A  the  stabilizing  and  destabilizing  terms  in 
Eq.  (1)  just  balance  each  other.  The  beam  has  no  free  energy  available 


for  driving  the  instability  when  y = y . , where  y . ■ 

J oi  i.crit  i,crit 

(1  - B2  . )-1^2.  The  roots  of  Eq.  (2)  can  be  easily  evaluated  in 
i.crit 

two  regimes:  (i)  B >8  . , and  (ii)  B >>  B . . In  the  first 

oi  i,cnt  oi  i.crit 

regime,  (i) , the  perpendicular  energy  of  the  beam  is  close  to  the  critical 


value  and  the  saturation  is  caused  by  the  depletion  of  free  energy.  In 
the  second  regime,  (ii),  the  perpendicular  energy  of  the  beam  is  well 
above  the  critical  value  and,  as  we  will  show  in  Section  VI,  the  satura- 
tion mechanism  for  the  instability  is  phase  trapping  of  the  particles 
by  the  excited  wave. 

Case  (1) 

Here,  8 ^ B . _ and  the  second  term  on  the  right  hand  side  of 

oi  i,crit 

Eq.  1 is  larger  but  comparable  to  the  first  term.  The  linear  frequency 
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shift  and  growth  rate  take  the  form: 

Aw  = Re(5wo).  = ((3/2)Znji  A2wq)  1/3  Awq,  (3a) 

T,  = Im(6w  ),  = 3~l/2(Aw  )2/3  (3Z  „/2)~2/3  ( (3Z  „/2)2  - A2w  )l/2 

1 o i q n*.  nx.  o 

= (Z  . Aw  - 2 A2w  /3)1/2  (3b) 

nil  o o 

The  growth  rate  for  this  case  is  proportional  to  the  fourth  root  of  the 
surface  charge  density. 

Case  (ii) 

Here,  8 >>  8 . and  the  second  term  on  the  righthand  side  of 

oi  i.crit 

Eq.  (2)  dominates.  The  linear  frequency  shift  and  growth  rate  of  the 
wave  are: 


Aw . . “ R (6w  ) = 

n e o ii 

(1/2)  <3Znt  A0>‘/S. 

(4a) 

T. . ” Im(6w  ) = 

11  O ii 

3 1 /2/2  (3  Z A2uj  ) 1/3. 

nx  o 

(4b) 

The  growth  rate  for  this  case  is  proportional  to  the  third  root  of  the 
surface  charge  density,  o . 

The  condition  that  the  nth  waveguide  mode  and  the  synchronous 

mode  intersect  on  the  dispersion  curve  limits  the  value  of  x . It  is 

n 

easy  to  show  that  x^  is  bounded  by  0 ^ ^ l.  Figure  3 depicts  the 

functions  and  as  functions  of  x^  for  l = 1,  2,  3 and  A.  Fig.  4 

shows  the  dispersion  diagram  for  the  cyclotron  maser  instability  when 

v » 0,  uj  = ck  * in  /y  and  S,  ■ n - 1.  Note  that  the  cyclotron  mode 
on  o n o o 

splits  into  two  branches  of  negative  and  positive  energy  polarity.  The 
negative  energy  cyclotron  mode  couples  to  the  positive  energy  TE  guide 
mode  resulting  in  an  instability.  The  growth  rate  maximizes  at  k^  - 0. 
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4 


If  the  cyclotron  frequency  were  greater  than  the  guide  cut-off  frequency. 


ft  / y > ck  , the  growth  rate  would  have  two  peaks  at  !k  ! >0,  sym- 

o ox  n 1 o 

metric  about  the  k =0  axis.  The  maximum  growth  rate  for  ft  /y  > ck 
o o o-t  n 

is  always  less  than  the  maximum  growth  rate  for  ft  /y  = ck  . This  can 

o ox  n 

be  seen  in  Fig.  5,  where  growth  rate  contours  are  shown  as  functions  of 


(ft  /y  )/(ck  ) and  (ck  /ck  ).  We  note  that,  in  general,  the  growth  rate 
o ox  n on  o.o 

is  maximum  when  the  wave  group  velocity,  v , equals  the  axial  beam 


velocitv,  v . This  becomes  apparent  if  we  transform  to  the  beam  frame 
on 

(v  * 0)  and  note  that  for  fixed  8 and  k , the  growth  rate  is  largest 


when  k = 0 (i.e.,  when  w = Jtft  /y  = ck  ) . Atk  *0,  the  group 
o o o o n o ’ 6 K 

velocity  vanishes  and,  therefore,  v = v * 0 in  this  frame.  Since 

g on 

both  v and  v transform  in  the  same  way,  the  growth  rate  maximizes 

g o II 

when  v = v in  all  frames, 
g on 
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Section  IV.  Nonlinear  Formalism 


In  this  section  the  basic  equations  governing  the  nonlinear  behavior 
of  the  electron  cyclotron  maser  instability  are  derived  using  a single 
wave  model.  The  validity  of  the  single  wave  approach  will  be  examined 
in  Section  V. 

In  our  procedure,  we  express  the  nonlinear  dynamics  of  a single 
wave  in  terms  of  an  ensemble  average  of  the  nonlinear  particle  orbits. 

The  particle  orbits  are  related  to  the  fields  through  the  relativistic 
Lorentz  force  equations.  The  resulting  wave  and  force  equations  are 
then  numerically  solved  to  obtain  the  self-consistent  behavior  of  the 
particles  and  the  field.  This  analysis  describes  the  nonlinear  evolu- 
tion of  the  field  amplitude  and  frequency  as  a function  of  time.  The 
beam  and  waveguide  configuration  shown  in  Fig.  2 is  chosen  as  our  basic 
model . 

We  have  noted  in  Section  III  that  the  linear  growth  rate  for  the 

excited  TE  mode  maximizes  when  the  axial  beam  velocity  and  the  wave 
on 

group  velocity  are  equal.  We  chose  to  analyze  this  situation  since  we 
are  interested  in  the  regime  of  maximum  wave  growth.  Furthermore,  for 
convenience,  we  transform  to  a frame  in  which  both  the  axial  beam 
velocity  and  wave  group  velocity  vanish.  Since  the  group  velocity  is 
now  zero  in  this  frame  (cutoff  frame),  we  note  that  the  axial  wave 
number  of  the  wave  is  also  zero.  The  electric  field,  in  the  cutoff 
frame,  of  the  vacuum  waveguide  mode  has  the  form 

Ey (x, t)  - - Eq (t)  cos  (w^t  + a (t) ) sin  (kn(x  ~ a)),  (5a) 
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where  k^  = 7rn/(2a)  and  n is  a positive  integer  such  that  E^(x,t)  satis- 
fies the  appropriate  boundary  conditions  at  x = ± a.  In  Eq.  (5a)  the 

frequency  uj  is  a constant,  while  the  amplitude,  E (t),  and  the  fre- 

o o 

quency  shift,  3oc(t)/3t,  are  weak  functions  of  time  (e.g.,  3 In  (E  (t))/3t, 

o 

3a(t)/3t  <<  uq) . The  amplitude  EQ(t)  can  be  expressed  as 


£ (t)  « G exp  f r ( t ’ )dt' , 
o o J 


where  c is  the  initial  field  amplitude  and  T(t)  is  the  time  dependent 
growth  rate.  We  note  that  in  the  linear  regime  of  the  instability  the 
frequency  shift  and  growth  rate  are  constant  and  equal  to  their  linear 
values,  3a(t)/3t  = AwT  and  F(t)  = T.  . 

L la 

Associated  with  the  electric  field  in  Eq.  (5a)  is  the  vector 


potential  A^(x,t)  given  by 


Ay (x,t)  = (c/w 


„>  I*1-1 


/ ol>  ) E sin  (w  t + a) 
o o o 


+ (E  /w  ) cos  (u  t + a)r  sin  (k  (x  - a)),  (6) 

oo  o ) n 

where  A (x,t)  is  valid  to  first  order  in  the  small  parameters  r(t)/uj 
y o 

and  i(t)/w  . The  time  evolution  of  E (t)  and  o(t)  is  determined  by  the 
oo 

particle  current  density,  Jy(x,t),  through  the  wave  equation  for  Ay(x,t): 


(32/3x2  - c“2  32/3t2)  A - - Arc-1  J . < 

y y 

Substituting  (6)  into  the  wave  equation  and  keeping  only  lowest  order 


terms  in  the  small  parameters  and  a/wQ  yields  the  relation: 
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•f(u)2  - c2k2  + (c2k^  + w2)( a/tu  ))  E sin  (ai  t + a)  • 
o n nooo  o 

- (c~k2  + (ii^)  (E  / u>  ) cos  ((i)  c + ct)>  sin  (k 

nooo  o I n 


(x  - a)) 


4iru > J (x,t). 
o y 


(8) 


The  current  density  for  a discrete  set  of  particles  is  given  by: 

N 

|e|o  ^ 

Jy(x,t)  — / , vy(4)i,t)  £(x  - x(f>i,t)),  (9) 


i=l 

where  is  the  surface  number  density  in  the  y-z  plane,  N is  the  number 
of  particles  on  a gyro  radius,  $ is  the  initial  velocity  space  angle 
of  the  iCtl  particle  measured  with  respect  to  the  x axis  and  x(#  ,t)  is 
the  x position  of  the  ith  particle.  In  the  limit  as  N + ", 

2lt 

N_1  2 "■  (2t)_1  f d*o 

o 

where  replaces  In  Eq.  (9)  the  velocity  component  v of  the  i 

particle  is  a function  of  only  ($  ,t)  since  we  assume  an  initially  cold 
particle  distribution  in  velocity.  Our  theory  can  include  an  initially 
thermal  particle  distribution,  which  would  add  an  average  over  the  ini- 
tial velocities  in  (9).  We  will  show  later  in  Eq.  (25)  that  for  an 

initial  energy  spread,  6y  <<  2y  Au)/(fi  /y  ),  the  particle  distribu- 

X ox  o ox 

tion  can  be  considered  monoenergetic  at  t * 0.  The  present  analysis 

is  valid  when  the  inequality  in  Eq.  (25)  is  satisfied.  This  does  not, 

however,  prevent  the  particles  from  thermalizing  as  the  instability 

develops.  The  wave  equation  can  be  separated  into  equations  for  E (t) 

o 

and  a(t).  By  multiplying  Eq.  (8)  by  sin  (k^(x  - a)),  integrating  across 
the  waveguide  from  x » -a  to  x » a and  operating  on  the  resulting 


) 
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equation  with 


t+2iT/(u>  -hi)  > 

/•  0 I sin  (b  t'  + a(t')) 

/ ’ ° ! 


) u \ 

| cos  (u)Qt ' + a (t ' ) ) J , 


we 


arrive  at  the  following  expressions,  which  are  valid  to  order 


r/u  , ot/oo  : 
o o 


“(t)  “ (U2  + c2k^) 


o n 


2 | e | a 

c2k2  - u2  + ) to  (w  + ct) 

no  /a  oo 


t+2n/(a)  +i) 


/ 


(Wq+oi)  | 

dt'  <v.(<t>„  t')  sin  (kn(x(<})o,t')  - a))  sin  (u^t'  +a(t'))>[  . 


y °> 


i 


f 


t-*-2Tr/ (to  +a) 


i “1 


It 


j*  E(t')  sin2  (u)ot'  +a(t'))  dt 


(10a) 


t+2ir/(u  +a) 


2 | e | 0q 

1 " u)2  + c^k2  a “o'  o 

o n 


(un  + a)  • / dt  <vy(*0.t  ) 


sin  (kn(x  (<|>o , t * ) - a))  cos  ^t’  + a(t')j>^>. 


(t+2Tf/(oj  -hi) 
o 


E(t')  cos2  (u  t'  + a(t'))  dt'  > 


(10b) 
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in  the  above  equations,  < > 


is  the  ensemble 


2 ji 

= (2n)“  1 J*  d:?o 
o 

average  over  the  initial  particle  phase.  We  note  that  the  temporal 

averages  in  Eqs.  (10a)  and  (10b)  are  performed  over  the  actual  wave 

period,  2 tt  / (co  + a).  Equations  (10a)  and  (10b)  describe  the  linear  as 
o 

well  as  nonlinear  evolution  of  the  wave  frequency  and  amplitude  in  terms 

of  the  particle  orbits.  The  only  restriction  placed  on  F(t)  and  a(t)  is 

that  they  vary  little  in  a wave  period,  2ir/((oQ  + a). 

The  righthand  side  of  Eqs.  (10a)  and  (10b)  contain  the  details  of 

the  particle  dynamics  which  are  related  to  the  fields  through  the 

relativistic  orbit  equations.  Defining  0 = 0 + i0  where  0 = v /c 

x y xx 

and  0 = v /c,  the  relativistic  orbit  equation  can  be  put  into  the  form: 


dS/dt 


i 

0 + - 


Bz(x,t) 


m c 
o 


S 


“^4  (1  - 0 (6  - 6*)/2)  Ey(x,t) 


where  y = (1  - 60*)-1^2,  E (x,t)  = - c-1  3A  /3t,  B (x,t)  = 

i y y z 

and  x is  given  by  dx/dt  - c(B  + 0*)/2. 

We  now  choose  the  solution  for  Eq.  (11)  to  be  of  the  form: 

S(4>  ,t)  = 0 (4>  ,t)  exp  (i 4>(4>  ,t)) 

O J.  o o’ 

where  0.  (4>  ,t)  = (02  + 02)  1/2  = (6B*)l/z  and 
■*•0  x y 

t 


«(<>  ,t) 
o 


■/ 


4; (♦  ,t')  dt'  + $ = tan 


-1 


(0  /6  ) 
y x 


3A  / 3x 

y 


(ID 


(12) 


is  the  velocity  space  angle  at  time  t.  Substituting  Eq.  (12)  into  (11) 
and  equating  real  and  imaginary  parts,  we  obtain  the  following  fully 
relativistic  equations  for  the  velocity  magnitude  and  velocity  space  angle 
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(13a) 


dg  — | e | E (x,t) 

~rr  = r^T” : sin  $($  ,t), 

dt  YJm  c o 

a o 


d_J 

dt 


|ej  E (x,'.) 
v 

Y tn  cB 

A O X 


| e , Bz(x,t) 

y m C 

± o 


cos  $($  ,t), 
o 


(13b) 


where  y =»  (1  - b2)-1/2  and  dx/dt  = c8  cos  $ (<j>  ,t).  Equations  (5),  (10) 

1 1 ± o 

and  (13)  form  a set  of  coupled  nonlinear  equations  which  describe  the 
evolution  of  the  electron  cyclotron  maser  instability  in  the  single  wave 
model . 

In  Appendix  A we  show  that  the  results  from  linear  theory  can  be 
recovered  from  these  equations.  The  orbit  equation  in  (13)  when  written 
on  the  slow  time  scale  yields  a constant  of  the  motion  which  is  dis- 
cussed in  Section  VI.  The  constant  is  useful  in  forming  a qualitative 
picture  of  the  trapping  and  saturation  process. 

The  amplitude  of  the  electric  EQ(t/  can  be  related  to  average 
decrease  in  particle  energy  through  the  energy  conservation  equation. 
Conservation  of  total  energy  within  the  waveguide  implies  that: 


jz  I (e,(x,t)  + e (x,t))  dx  » et. 


(14) 


where  ef(x,t)  ■ (E2(x,t)  + B2(x,t))/8ir  is  the  field  energy  density, 
t y z 

£ (x,t)  - o <(y(<t  t)  - 1)  6 (x  - x (<})  ,t))>  m c2  is  the  particle 
poo,  o o 

energy  density  and  e.  - (o  /2a)  (y  - 1)  m c2  is  the  average  total 

to  o o 

energy  density.  Substituting  E - - c-1  3A  /3t  and  B - 3A  /3x  into 

y y z y 
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(14),  where  A^  is  given  by  Eq.  (6),  and  performing  the  spatial  average 
as  well  as  a temporal  average,  leads  to  the  following  expression  for 


E (t): 
o 


■1/2 


2 m c ck  \z 

E (t)  = “ITT  “h  1 + — <1  " 2 «/“  ) 

o |e|  b \^0/  0 


.1/2 

( y - <y  ($  ,t)  >)  , 

\ Oi.  J.  o / 


(15) 


where  a>.  = 4t  ' e j 2 a / (m  a)  is  a modified  plasma  frequency  and  y is 
b 1 o o oj. 

the  initial  gamma  of  the  particles.  In  deriving  Eq.  (.15),  terms  of  order 


higher  than  T/u  and  a/w  have  been  neglected.  Since  ck  /u  1 and 
o o no 

a/oj  <<  1,  the  expression  for  E (t)  can  be  simplified  to: 
o o 


/§  m c 

E (t)  * | -r°-  «. 

o e b 


v ' < Wc)  > 


i/: 


(16) 


A simple  analytic  expression  for  < y ^ (j  , t) > at  saturation  which  will 
permit  us  to  determine  the  maximum  fr’eld  strength  through  Eq.  (16),  is 
obtained  in  Section  VI. 

Since  our  system  of  fields  and  particles  is  independent  of  the 

spatial  variable  y,  we  note  that  the  y component  of  canonical  momentum 

P » y m v - c-1  Iel  A (x,t)  - 0 m x is  a constant  of  the  motion  for 
y a o y 1 1 y o o 

each  particle.  The  conservation  of  total  energy  as  well  as  the  y 
component  of  canonical  angular  momentum  is  monitored  throughout  the 
numerical  simulations  of  Eqs.  (10)  and  (13)  to  ensure  consistency. 

The  nonlinear  model  developed  in  this  section  is  directly  applicable 
to  the  experimental  situation  in  which  a wave  is  spatially  amplified  in 
the  streaming  direction  of  a gyrating  beam.  The  saturation  levels  ob- 
tained with  the  present  model  are  directly  related  to  the  saturation 
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levels  in  an  amplifier  when  the  axial  beam  velocity  equals  the  wave 
group  velocity.  As  pointed  out  in  Section  III,  this  situation  corres- 
ponds to  a grazing  intersection  of  the  waveguide  mode  with  the  beam 
cyclotron  mode  when  viewed  on  the  dispersion  diagram.  Furthermore, 
such  a grazing  intersection  leads  to  an  absolute  maximum  in  the  temporal 
growth  rate.  In  a steady  state  amplifier  the  input  wave  amplitude 
initially  grows  exponentially,  enters  the  nonlinear  regime  and  saturates. 
If  we  follow  a group  of  particles  contained  in  a thin  cross-sectional 
slab  of  the  beam,  we  note  that  the  net  flow  of  total  energy  flux  into 
this  slab  is  zero  since  we  only  consider  the  case  where  the  axial  beam 
velocity  and  wave  group  velocity  are  equal.  Therefore,  in  the  reference 
frame  of  the  particles,  the  field  amplitude  evolves  in  time.  This  is 
precisely  the  situation  described  by  our  model.  With  the  appropriate 
transformations  (see  Appendix  B) , the  saturation  levels  of  the  spatially 
amplified  wave  can  be  obtained  with  the  present  temporal  model. 
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Section  V.  Nature  of  Approximations 


A number  of  approximations  have  been  made  in  develop!:.^  > r 

linear  model.  Among  these  are  i)  the  assumption  of  a single 

wave,  ii)  the  assumption  of  a single  excited  vacuum  mode  in  th. 

guide  and,  finally,  iii)  the  initially  monoenergetic  beam  appr  ximu : . 

In  this  section  the  regimes  of  validity  of  these  approximations  are 

examined.  We  find  that  our  model  is  valid  for  a wide  range  of  parameters. 

To  verify  the  validity  of  the  single  wave  approximation  the  growth 

rate  and  field  energy  spectrum  are  examined.  Using  a Taylor  series 

expansion  of  the  linear  growth  rate  about  = 0 , r(kQ)  ^ 1(0)  + 

(1/2)  (32T(0)/3k2)  k",  we  find  that  the  half  width  of  the  growth  rate 
o o 

spectrum  is  Ak  = 2 (T (0) / | 32T (0) /3k2 | ) 1 ^ . Since  the  field  energy  is 

proportional  to  exp  (2T(0)t),  the  half  width  of  the  energy  spectrum 

after  the  field  has  e-folded  N times,  t = N/r(0)  is 

6k  = 2 (In  2/N)1/2  (T (0)/ I 32T (0) /3k2 | ) 1 ! 1 = N-1/2  (In  2)l/2  Ak.  (17) 

From  the  linear  dispersion  relation,  Eq.  (1),  we  find  that,  for  the  two 

cases  discussed  in  Section  II,  the  half  width  of  the  linear  growth  rate 

spectrum  for  Case  (i)  (i.e. , 8 - 8 . „ and  w = ck  = in  / y ) 

ox  x,crit  o n o ox 


Ak  % 2/6  (I 


/u.  \‘/2 

‘ <0>/c)  te) 


(13a) 


and  for  Case  (ii)  (i.e.,  8 >:>  8 . and  to  = ck  = m /y  ), 

ox  x,crit  o n o ' ox 


Ak . . 5:  to  /c  , 
ii  o 


(13b) 


In  Eq.  (18a),  [^(0)  is  the  linear  growth  rate  at  k^  = 0 as  given  by  Eq. 
(3b).  The  half  width  of  the  energy  spectrum  after  the  field  ampli-.ude 
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has  e-folded  N times  is  for  Case  (i) 

, tjj  \l/2 

5k.  * 2.6  (In  2)1/2  N-1/2  (T . (0)/Au>  )(  — -)  <u  /c  (19a) 

i 1 ° \ “c  ) o 

and  for  Case  (ii) 

6k  . = (In  2)1/2  N_1/2  u>  /c  (19b) 

ii  o 

The  frequency  spread  corresponding  to  the  wavenumber  spread  in  Eqs.  (19a, 

b)  can  be  estimated  from  the  vacuum  dispersion  relation  ui  = (w  2 + 

o 

, ,1/2 

c^k^-)  . The  frequency  spread  is  given  approximately  by  6tu  =* 

(c6k/^Q) . From  the  expressions  for  the  half  widths  of  the  energy 

spectrum  in  Eqs.  (19a,  b)  and  the  frequency  spread,  we  see  that  the 

spectrum  can  be  represented  by  a single  wave  if 

N-!/2  (r  (0)Mu>  ) (Aw  /u  W2  « 1 
1 O O o 

for  Case  (i)  and  N_1/^2  <<  1 for  Case  (ii).  Typically  in  the  examples 
we  shall  consider  the  number  of  e-folds  before  the  wave  saturates  is 
N * 15. 

In  the  nonlinear  analysis,  the  general  form  of  the  excited  wave  was 
that  of  a single  vacuum  waveguide  mode.  In  the  presence  of  a tenuous 
electron  beam,  this  vacuum  mode  is  expected  to  be  a good  representation 
of  the  actual  field  in  the  loaded  guide.  If  the  particle  density  is 
sufficiently  large,  the  field  induced  by  the  beam  current  density  will 
couple  to  more  than  one  vacuum  waveguide  mode.  Consequently,  the  pro- 
file of  the  actual  field  will  differ  from  the  form,  sin(k  (x-a)),  used 

n 

in  the  analysis.  We  now  show  that,  in  the  tenuous  beam  limit,  the 
correction  to  the  vacuum  waveguide  field  is  of  order  Au/uq  and  can, 
therefore,  be  neglected. 
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To  obtain  a rough  estimate  for  the  correction  to  the  vacuum  field 


we  assume  that  the  particle  current  density  at  the  frequency  tu  is 
produced  solely  by  the  vacuum  mode  given  in  Eq.  (5).  To  simplify  the 
calculation  we  take  x <<  a.  The  electric  field  in  the  loaded  waveguide 
is  chosen  to  have  the  form 


E (x,t)  = - E (t)  cos  (uj  t + a(t))  S a sin  (k  (x  - a)),  (20 

y o o ' m m 

m=l 


where  E (t)  , u)  , a(t)  and  k have  been  defined  following  Eq.  (5)  and 
o o m 


a are  the  constant  relative  amplitude  coefficients  of  the  electric 
m 


th 


field.  If  corrections  to  the  n vacuum  mode  are  desired,  we  assume 

th 


that  the  current  density  is  produced  by  the  n vacuum  mode  given  in 
Eq.  (5).  Substituting  the  vector  potential  associated  with  Eq.  (20) 
into  the  wave  equation,  (7),  we  find: 


m=l 


- c2k-  + (c‘k‘  + oj^)  d/u) 
o m moo 


E (t)  sin  (u  t + a) 
o o 


- (c2k/  + in 2)  (E  /a)  ) cos  (to  t + a) 
m o o o o 


> a sin  (k  (x  - a)) 
m m 


= - A 7TW  J (x , t)  . 
o y 


(2i: 


The  current  density  in  Eq.  (21)  can  now  be  written  as: 


J (x,t)  “ J (t)  sin  k (x  - a) , 

y m m 

ra*l 


(22) 


where 


f , (nH-l)/2 

I -(  | e | o / (2ira)  ) (-  1)  <v  ($  ,t)>,  for  m odd 

j (t)  - < 0 y 0 

ra  ) 0 


otherwise 


when  k x($  ,t)  <<  ~/2.  Combining  Eq.  (21)  and  (22)  and  neglecting 
m o 

higher  order  terms,  we  find 


m 4 y D~ 

04- 


'(1+7  2) 

_ °X 

^l 


- e2  (0,2  - c:k2) 1 

oi  7 D *•  / m 

oi 


(23) 


where  m ^ n and  the  linear  expression  for  v^^.t)  given  by  Eq.  (A4) 
was  used.  Comparing  the  right  hand  side  of  (23)  with  (A5a)  shows  that: 


((oj“-  - c^k-)  + (u~  + c2k“)  Aw/w) 
a ^ o n o n 

m 


CO 


2 - c2k“- 


m 


which  reduces  to 


2 (Aco/co  ) 

a ^ i <<; 
m 1 - (m/n)z 


(24) 


since  5:  ck^.  The  coupling  to  other  vacuum  waveguide  mode  is  of  order 

Aco/co  and  hence  our  original  choice  for  the  electric  field  given  in  (5) 
o 

is  indeed  reasonable. 

In  this  paper  we  consider  the  beam  to  be  initially  monoenergetic  in 
the  beam  frame.  The  range  of  validity  of  this  assumption  requires  some 
examination.  Since  reference  will  be  made  to  quantities  in  both  the 
laboratory  and  beam  frames,  we  shall  denote  the  former  with  unprimed 
variables  and  the  latter  with  primes.  It  is  clear  that  the  cold  beam 
approximation  will  be  valid  if  the  total  initial  spread  in  yJ.  satisfies 
<5 <<  2 (y^  ~ <Y^>S>»  where  < Y^  >s  is  the  average  value  of  gamma 

at  saturation.  In  terms  of  efficiency,  this  condition  becomes 


5y’/y'  <<  2 n' 

i oi 


(y1 

oi 


D/y’  , 

oi 


i a 


(25) 


where  n'  is  the  beam  frame  efficiency  defined  as  p'  = (y^a  “ 

(y  - 1).  Since  the  unstable  wave  has  an  infinite  parallel  wavelength 

ox 

in  the  beam  frame,  any  spread  in  v,J  which  does  not  violate  the  above 
inequality  is  unimportant. 

To  define  the  beam  qualities  necessary  for  the  validity  of  our  cold 
beam  approximation,  it  is  useful  to  express  Eq.  (25)  in  terms  of 
laboratory  frame  quantities.  The  righthand  side  of  Eq.  (25)  can  be 
written  in  terms  of  laboratory  frame  quantities  by  using  the  trans- 
formation rule  for  efficiency  derived  in  Appendix  B,  i.e., 

n = n'  (y * - 1)  Y /(y'  (y  - 1))  where  n is  the  power  efficiency  in 

o-*-  o ox  o 

the  laboratory  frame.  Using  the  results  of  Appendix  B,  the  inequality 
in  Eq.  (26)  becomes 


2 (Y„ 


1) 


<< 


n. 


(26) 


where  y = (1  - v2  /c^  - v2  /c2) 
o o II  ox 

A spread  in  the  beam  frame  gamma,  <5y^>  can  arise  from  (a)  a finite 
beam  temperature  in  the  laboratory  frame  or,  in  the  case  of  a mono- 
energetic  beam  in  the  laboratory  frame,  (b)  a spread  in  the  velocity 
space  angle  of  the  particles. 

In  case  (a)  only  a temperature  spread  in  the  perpendicular  velocity 

is  considered,  so  that  v * v + Av  where  Av  <<  v is  the  thermal 

x ox  x x ox 

component  of  velocity.  Only  perpendicular  temperature  is  considered 

since  a parallel  temperature  in  the  laboratory  frame,  AvM  , contributes 

a higher  order  correction  to  5y'  than  Av  if  Av„  and  Av  are  of  the 

x x " x 

same  order.  Since  the  perpendicular  momentum  is  frame  invariant. 
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yv  = y'v* , we  find  that  the  thermal  spread  in  y takes  the  form: 
x xx 

Ay  = ± X Ay'  (2:) 

Oil  -1- 

-1/2 

where  y = (1  - vz  /c2)  , v is  the  axial  laboratory  frame 

on  c n on  J 

velocity  and  Ay!  = 2(v  /c)  Y3  Av  /c  is  the  thermal  spread  in  y due  to 

Av'  * y Av  . Substituting  Eq.  (27)  into  (2b)  gives  the  following 

•l  on  ■*. 

requirement  for  neglecting  thermal  effects  in  vA: 

Ay/(yq  - 1)  « 2 n.  (23) 

In  the  above  inequality  the  efficiency,  p , is  calculated  using  the  cold 
beam  approximation. 

In  case  (b)  the  particles  have  the  same  energy,  however,  a spread 

in  the  velocity  components  exists.  Consequently,  the  particle  velocity 

components  in  the  laboratory  frame  can  be  written  as 

v = v - v A8,  (29a) 

r or  on 

v = v + v A6 , (29b) 

11  on  or 

where  | A0 | <<  1 is  the  spread  of  the  velocity  space  angle,  i.e.,  the 

angle  between  the  total  velocity  vectors  of  the  particles.  To  first 

order  in  A8 , we  see  that  v2  + v2  = v2  + v2  , so  that  the  system  of 

particles  is  monoenergetic  in  the  laboratory  frame.  The  transformation 

rule  for  y from  the  laboratory  to  the  beam  frame,  y'  - y y 

on 

(1  - v v /c2),  shows  that  the  particles  are  not  monoenergetic  in 
on  n 

the  beam  frame.  In  the  beam  frame  we  find  that: 


C + AY1’ 


(30) 
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where  y’  = y /y  and  Ay'  - y y (v  v /cMAG.  Substituting 
o-i-  o Oil  a on  o cii  oj. 

into  Eq.  (26)  gives  a condition  on  the  magnitude  of  the  maximum 
angular  spread. 


I aq  | 


max 


<< 


2 (^Q  - n 

y y*-  v V Tc~  * 

o Oil  Oil  Ol 


(31) 


allowed  in  the  cold  beam  approximation.  We  will  return  to  conditions 
(28)  and  (31)  in  the  next  section. 

One  of  the  causes  of  energy  spread  in  an  unneutralized  beam  is 
the  self  electrostatic  field.  For  the  planar  beam  that  we  have  con- 
sidered, the  total  spread  in  y ^ can  be  shown  to  be: 


5y 9 * ( »')2  Bfa<(0  /y\)/ck  ) y9 

x,self  oj.  o oj.  n oi 


(32) 


where  u*  * u ' / ( vy9  u9 ) and  u9  * ft  /y9  • Considering  Eq.  (25),  this 
b Ox  o o o 01 


imposes  an  upper  limit  on  p1,  and  hence  the  beam  current  in  the  labora- 


Section  VI.  Saturation  Mechanisms 


In  this  section  some  qualitative  arguments  concerning  the  saturation 
of  the  electron  cyclotron  maser  instability  are  made.  There  are  two 
mechanisms  which  are  responsible  for  saturation  of  the  unstable  wave  in 
the  maser  instability.  They  are  i)  depletion  of  the  rotational  free 
energy  of  the  electrons  and  ii)  phase  trapping  of  the  gyrating  particles 
in  the  wave.  Which  of  the  two  mechanisms  is  actually  responsible  for 
saturation  in  a particular  situation  depends  on  the  initial  beam 
parameters . 

Case  (i)  Free  Energy  Depletion 

As  mentioned  in  Section  III,  the  linear  dispersion  relation, 

Eq.  (1),  gives  rise  to  a threshold  condition  for  instability  given  by 
Boi  > 6a  cr^t«  That  is,  the  maximum  free  energy  per  particle  available 
to  the  wave  is 


e , * (y  - y . ) m c2 , 

free  ox  x,crit  o 


(33) 


where  y^  crit  * (1  ~ crit^~  ^ the  ^eam  Particles  were  to  lose 
all  of  the  free  energy  given  by  (32),  the  energy  conversion  efficiency 
would  be 


(y 


ox  'x.crit 


(>ox  - 


(34) 


where  yQA  ^ y ^ crit-  This,  of  course,  is  a rough  approximation,  since 
as  the  instability  develops  the  particles  warm  up,  while  the  threshold 
condition  was  obtained  for  a cold  beam.  Furthermore,  as  will  be  seen, 
competing  processes  such  as  electron  trapping  may  take  place  before  the 


beam  loses  its  free  energy.  However,  if  y is  slightly  greater  than 


Ya  cr^t»  the  expression  in  Eq.  (33)  will  represent  a good  approximation 

to  the  actual  conversion  efficiency.  Since,  if  v ^ y . , the 

oi  x,crit 

particles  lose  all  their  free  energy  before  thermalizing  or  trapping 

takes  place  and  the  wave  will  be  linearly  stabilized.  If,  however, 

Y >>  y . , the  particles  will  phase  trap  in  the  wave  and  satura- 

oi  x,crit 

tion  will  occur  before  all  the  free  energy  is  used  up. 

Case  (ii)  Phase  Trapping 

The  condition  for  wave  saturation  due  to  phase  trapping  can  be 

viewed  in  a number  of  ways.  First,  we  have  noted  in  Sections  II  and 

III  that  the  frequency  of  the  wave  must  be  slightly  greater  than  the 

relativistic  electron  cyclotron  frequency,  /y  » for  an  instability 

to  develop.  Therefore,  initially  we  have 

w - fl  /y  * Aw  ~ 0,  (35) 

o ox 

where  Aw  is  the  frequency  shift  which  can  be  obtained  from  the  linear 
dispersion  relation.  Now,  as  the  instability  develops  the  average  gamma 
of  the  particles,  <YA>  . decreases  until  w - Oq/  <yJL>  = - Aw.  At 
saturation  <y^>  is  minimum  and: 

w - Cl  /<Y  > - “ (36) 

o a s 

where  <ya>s  is  the  average  value  of  y^  at  saturation.  Comparing  (35) 

with  (36),  we  find  that  at  saturation  the  average  change  in  the  gamma  of 

the  particles,  <Ay  > “ Y - <y  > , is  given  by: 
x s ox  xs 

<Ay,  > “ 2y  Aw/w.  (37) 

x s ox 

The  efficiency  of  conversion  of  kinetic  energy  to  field  energy,  n,  when 
electron  phase  trapping  is  responsible  for  saturation,  is: 

n « 2 (Aw/w)  Y / (y  - 1).  (38) 

ox  ox 
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It  should  be  noted  that  the  efficiency  in  the  beam  frame  cannot  be  made 


arbitrarily  close  to  100  percent  simply  by  increasing  the  frequency  shift, 

■Aw.  The  growth  rate  of  the  field  is  a function  of  Aw  and  is  nonvanishing 

for  a limited  range  of  Aw.  The  relations  in  (37)  and  (38)  are  valid 

only  if  Aw  lies  within  the  nonvanishing  part  of  the  growth  rate  spectrum. 

To  find  the  actual  amplitude  of  the  field  at  saturation,  E , we  use 

o,s 

the  expression  for  conservation  of  total  energy  given  in  Eq.  (16).  The 
maximum  field  amplitude  at  saturation  is  easily  shown  to  be: 


m c 

E =4  -pn-  u>  (y  Aw/w  )1/,2» 
o ,s  ; e | b ox  o 


(39) 


where  it  has  been  assumed  that  a/w  <<  1 and  ck  /w  = 1. 

o no 

We  now  examine  the  nonlinear  behavior  of  the  particles  with  the 
aid  of  a constant  of  the  motion  derived  from  the  orbit  equations. 

The  orbit  equations  given  in  Eq . (13)  can  be  written,  in  the 
small  Larmor  radius  approximation,  x <<  a,  as: 


du 

i 

dt 

d£ 

dt 


- I e | E 

» a cos  ((w  + Aw)t)  sin  (p(t)  + $ ) 

m c o o 

o 

| e | E 

" — — cos  ((w  + Aw)t)  cos  (<p( t)  + p ) 

y . u m c o o 

i i o 


(40a) 

(40b) 


where  u±  * yA6A  is  proportional  to  the  transverse  particle  momentum, 

and  y ■ (1  + u2)1/2.  The  field  amplitude,  E , as  well  as  the  fre- 
11  o 

quencv,  w + Aw,  are  considered  fixed  and  w * fi  /y  is  the  initial 
o o o 01 

electron  cyclotron  frequency.  The  particle  momentum,  u m c,  and  time 

i o 

rate  of  change  of  the  velocity  space  angle,  dp/dt,  consist  of  fast  and 
slow  time  scale  variations. 
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The  fast  time  scale  variation  is  associated  with  the  wave  fre- 


quency or  electron  cyclotron  frequency,  while  the  slow  time  scale  is 
associated  with  the  frequency  shift,  Auk  The  nonlinear  behavior  is 
governed  predominantly  by  the  slowly  varying  parts  of  Uj_  and  d$/dt. 
Since  (wq  + Aw)  - d$/dt  changes  slowly  the  orbit  equations  in  (40) 


written  on  the  slow  time  scale  become 
du,  I e I E 


dt 


2 m c 
o 


sin  (X) 


(41a) 


,,  n e E 

% - (w  + Aw)  - -2-  4-  ' ° 

dt  o Y,  2 m c u. 


cos  (A) 


(41b) 


where  X = (wq  + Aw)t  - $(0  ~ $ • These  slow  time  scale  orbit  equations 
have  the  following  constant  of  motion: 


C 


c ° Ux(t)  cos  (X(t))  + (wQ  + Aw)  (yj.(t) 


2 

Y w \ 

OX  Q 1 

w + Aw/ 
o 


(42) 


Each  particle  moves  in  phase  space,  (uitX),  on  C = constant  curves. 

These  C * constant  curves  describe  the  particle  orbits  when  a constant 

frequency  and  constant  amplitude  field  is  instantaneously  switched  on. 

This  situation  can  be  realized  by  injecting  particles  into  a field 

filled  cavity.  The  particle  trajectories  in  phase  space  lie  on  the 

constant  C curves  depicted  in  Fig.  6.  The  initial  conditions  for  the 

example  shown  in  the  figure  are  Y„,  - 1.05,  p * 0.05,  Aw/w  - 0.018  and 

I e ! E / (y  w me)*  0.0025.  In  this  illustration  the  frequency  shift 

' 1 o cx  o o * * 

is  consistent  with  linear  theory  and  the  field  amplitude,  E , corres- 

o 

ponds  to  roughly  a quarter  of  the  value  at  saturation. 
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In  Fig.  6 the  particles  are  initially  uniformly  distributed 

between  0 ^ \ £ 2n  with  u^  = u . The  trajectory  followed  by  a 

particle  is  determined  by  its  particular  value  of  C which  is  different 

for  each  particle.  The  particles  have  access  to  a limited  region  of 

phase  space,  the  boundaries  of  which  are  shown  in  Fig.  6.  Particles 

which  initially  lie  on  closed  C * constant  curves  are  considered 

trapped  even  though  ux  as  a function  of  A is  single  valued.  However,  as 

the  particles  travel  along  their  C = constant  curve,  Uj_  will  eventually 

become  a multivalued  function  of  A.  Similarly,  for  a growing  field 

amplitude,  particles  can  be  considered  trapped  if  they  lie  on  closed 

C = constant  curves,  long  before  ux  becomes  a multivalued  function  of  A. 

The  character  of  the  accessible  C = constant  curves  is  determined,  among 

other  things,  by  the  value  chosen  for  E . For  sufficiently  small  values 

o 

of  Eq  the  particles  will  not  lie  on  closed  C = constant  curves  and, 

hence,  are  not  trapped.  The  value  of  Eq  for  which  the  first  closed 

C * constant  curve  intersects  the  u = u line  determines  the  field 

X ox 

amplitude  necessary  for  the  onset  of  trapping.  The  separatix  separating 

the  closed  and  open  orbits  first  intersects  the  line  u = u at  A = it. 

x ox 

As  larger  values  of  Eq  are  chosen,  the  inaccessible  interior  region  of 

phase  space  will  contract  and  approach  the  point  u = u and  A = . 

x ox 

At  the  value  of  E^  for  which  the  inaccessible  interior  region  vanishes, 
the  average  particle  energy  ceases  to  decrease.  The  value  of  electric 
field  amplitude  necessary  for  the  disappearance  of  the  interior  region 
is 


E = 
o 


2 


m cm 
o 


o 


I 


e i 


u 

ox 


Ato 

10 

o 


(43) 
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In  Figs.  7 and  8 the  actual  positions  of  the  particles  in  phase 
space  are  shown  at  various  times.  The  parameters  used  for  these  figures 
are  the  same  as  those  used  in  Fig.  6.  Superimposed  on  the  curves  in 
Figs.  7 and  8 are  the  boundaries  of  the  accessible  region  of  the  C = 
constant  curves.  At  t = 0,  when  the  electric  field  is  switched  on,  the 
particles  are  uniformly  distributed  along  the  u = u0  line  between 
0 £ X 2tt.  In  the  constant  amplitude  electric  field  the  particles 
proceed  to  bunch  and  rotate  about  the  interior  boundary  as  shown  in  the 
figures.  At  t * t^  the  particles  begin  to  overtake  each  other  and  bunch 
along  the  righthand  side  of  the  interior  boundary.  At  a still  later 
time,  t = t^,  the  particles  rotate  and  bunch  along  the  bottom  of  the 
interior  region.  At  this  point  the  configuration  corresponds  to  the 
saturation  state  since  the  average  particle  energy  is  at  a minimum.  The 
average  kinetic  energy  continues  to  oscillate  as  the  particles  rotate 
about  the  interior  boundary  in  a clockwise  direction.  The  distribution 
functions,  fe(uA),  corresponding  to  the  phase  space  trajectories  in 
Figs.  7 and  8 are  shown  in  Fig.  9.  As  can  be  seen  from  Figs.  8 and  9, 
not  all  the  particles  are  located  at  the  bottom  of  the  interior  boundary 
at  saturation.  Consequently,  the  average  of  u^  at  saturation  is  not 
exactly  given  by  uA  at  the  bottom  of  the  interior  region. 

For  sufficiently  small  values  of  EQ  the  accessible  C = constant 

curves  are  open  and  particles  are  not  trapped.  As  E is  increased, 

o 

the  first  particle  to  become  trapped  is  always  the  one  initially  located 
at  (uqa  it).  Trapping  begins  when  the  particle  at  (u^.tt)  has  a closed 
orbit  with  turning  points  at  A = 2tt  and  0.  At  these  turning  points, 
dA/dt  vanishes.  Since  the  constant  C is  invariant  along  this  closed 
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orbit  and  dX/dt  vanishes  at  X = 2tt  and  0,  we  can  obtain  the  threshold 

value  of  E for  trapping  from  Eqs.  (41b)  and  (42).  The  approximate 
o 

value  of  the  field  amplitude  for  the  onset  of  trapping  is  given  by: 


Using  the  constant  of  the  motion  in  Eq.  (42),  the  minimum  value  of  ut 
reached  by  the  first  trapped  particle  occurs  at  X = it  and  is  approxi- 
mately: 

u = u - ( —\.  <«) 
i.rnin  ox  uQi  \ uQJ 


The  corresponding  minimum  gamma  for  this  particle  is  Y ox 

2 v Aw  /w  . Consequently,  the  first  trapped  particle  undergoes  a 

Ojl  0 0 

change  in  gamma  at  1 = it  given  by: 


ay 

x 


Y Aw  /w  . 
'oi  oo 


(46) 


This  expression  for  the  maximum  change  in  gamma  allows  us  to 
estimate  the  field  amplitude  at  saturation  as  was  done  in  Eq.  (39).  The 
change  in  gamma  given  in  Eq.  (46)  is  identical  to  the  value  found 
in  Eq.  (37),  where  a different  line  of  reasoning  was  used. 
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Section  VII.  Results 


The  nonlinear  coupled  equations  (10a, b)  and  (L3a,b)  are  solved 
numerically  for  a wide  range  of  parameters.  In  each  case  the  distri- 
bution function  was  represented  by  approximately  40  to  100  test 
particles.  The  actual  number  of  particles  used  was  determined  by 
performing  runs  with  larger  and  larger  numbers  of  particles  until  no 
discernible  difference  appeared  in  the  results.  The  conservation  of 
both  total  system  energy  and  canonical  y-momentum  was  monitored.  The 
integration  timesteps  were  adjusted  to  maintain  the  conservation  of 
both  constants  to  within  0.5  percent  over  the  duration  of  a run.  As 
diagnostic  aids,  plots  of  the  test  ^article  velocity  distribution  func- 
tion, velocity  space  and  phase  space  were  made  at  various  times  during 
the  runs. 

A run  was  initiated  by  uniformly  distributing  the  particles  in  the 
velocity  space  angle,  0 ^ $ » 2 tt  and  assigning  them  the  same  perpen- 

dicular velocity,  v^  = v^.  A small  amplitude  electric  field  was 
introduced  as  a perturbation  and  allowed  to  grow  self  consistently. 

From  the  results  of  the  single  wave  simulation  runs  we  have  made  com- 
posite graphs  of  the  transfer  efficiency  from  particle  kinetic  energy 
to  wave  electromagnetic  energy  at  saturation  as  a function  of  several 
parameters.  A comparison  between  these  results  and  the  analytic  pre- 
dictions in  Section  VI  shows  good  qualitative  agreement. 

We  will  now  discuss  the  phase  space  dynamics  for  the  two  saturation 
mechanisms  discussed  in  Section  VI.  The  evolution  of  the  particles  in 
phase  space  and  their  distribution  functions  is  shown  when  saturation 


3c 


occurs  by  (i)  free  energy  depletion  and  (ii)  phase  trapping.  In  the 

free  energy  depletion  case  the  initial  parameters  are  yQ  = 1.012  ~ 

\ = 1.0108,  a = 0.05,  where  u is  a dimensionless  parameter  re- 

i.crit 

lated  to  the  beam  density  through  u = u,  /(/y  u ) , w = ck  = if.  / y 

boxoo  n o Ox 

and  l = n = 1.  Figures  10(a)  and  10(b)  show  the  particle  phase  space  in 
the  course  of  the  instability  when  the  growth  is  still  linear.  At  satura- 
tion, shown  in  Fig.  10(c),  the  particles  indicate  no  phase  trapping.  The 
average  particle  f A at  saturation  is  within  a few  percent  of  y ^ cr^t- 
However,  the  individual  particles  energies  have  a large  variation,  as 
can  be  seen  from  the  distribution  functions  depicted  in  Figs.  11(b)  and 
11(c).  The  predicted  efficiency  according  to  Eq.  (34)  is 

n = (,  -v  . )/(y  - 1)  ~ 10  percent.  This  is  in  excellent  agree- 

ci  x ,crit  ox 

ment  with  the  numerically  evaluated  efficiency  of  10  percent.  To 
illustrate  case  (ii)^  where  phase  space  trapping  is  responsible  for  satu- 
ration, we  choose  the  parameters  y = 1.2,  u = 0.05.  (Note  that 

Y . = 1.01  for  this  case  as  in  case  (i)).  The  particle  dynamics  in 

'x,crit 

velocity  phase  space  and  the  velocity  distribution  function  are  shown 
at  various  times  T during  the  evolution  of  the  run  in  Figs.  12  and  13. 

In  Fig.  12(a)  the  particles  are  shown  at  T =250  when  they  have  begun  to 
display  some  bunching  but  are  relatively  close  to  their  initial  phase  space 
positions.  The  corresponding  distribution  function  in  Fig.  13(a)  displays 
little  spread  in  velocity.  When  the  nonlinear  effects  manifest  them- 
selves by  a change  in  the  growth  rate,  the  velocity  phase  space  bunching 
becomes  distinctly  visible  and  the  velocity  distribution  begins  to 
spread  as  shown  in  Figs.  12(b)  and  13(b).  In  Figs.  12(c)  and  13(c),  we 
show  the  velocity  phase  space  and  distribution  function  at  saturation. 
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The  phase  space  plots  in  Figs.  12(b)  and  12(c)  show  the  trapped  electrons 
circulating  and  forming  a tongue  configuration,  similar  to  the  one  in 
Figs.  7 and  8,  where  a constant  electric  field  was  present.  The  satu- 
ration mechanism  for  this  case  is  clearly  phase  trapping  of  the  beam 
electrons  in  the  wave.  Because  not  all  the  electrons  are  trapped,  an 
average  over  the  distribution  function  is  necessary  to  determine  <YX>S- 
This  necessitates  the  introduction  of  a numerical  factor,  f,  in  the 
simple  expression  corresponding  to  Eq.  (37) 

<Ay  > /y  = 2f  Aw/w 

2L  S O*  O 

The  factor  f is  a number  on  the  order  of  unity  and  depends  on  the  frac- 
tion of  deeply  trapped  particles,  f ^ 1. 

As  an  example  of  a typical  run,  we  shall  follow  the  evolution  of 
field  and  particle  quantities  as  a function  of  time.  The  initial  con- 
ditions chosen  for  this  example  are:  y = 0.05,  y * 1.05,  w = ck  = 

oi  o n 

fi  /y  and  n = 1,  and  y * w. /(/y  u>  ) . In  Fig.  14,  the  nonlinear  fre- 
o oj.  b oi  o 

quency  shift  and  growth  rate  are  shown  as  functions  of  the  normalized 

time  parameter  x = UJ0t*  Initially  for  T ^ 150  when  the  field  amplitude 

is  small,  we  see  a transient,  during  which  the  perturbing  field  comes 

into  equilibrium  with  the  particles.  That  is,  initially  the  presence 

of  the  small  amplitude  electric  field  is  inconsistent  with  the  initial 

conditions  of  the  beam  particles.  During  the  early  transient  the  field 

and  particles  adjust  themselves  to  self  consistently  satisfy  Eqs.  (10a, b) 

and  (13a, b).  The  transient  regime  lasts  for  a time  roughly  given  by 

T ~ w /T  . The  linear  regime,  150  ~ T * 350,  follows  the  initial 
tran  o L 6 
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transient.  In  this  regime  the  wave  grows  exponentially  and  both  the 
growth  rate  and  frequency  shift  assume  the  values  predicted  by  linear 
theory.  As  will  be  apparent  from  one  of  our  later  figures,  the  growth 
rates  and  frequency  shifts  in  the  linear  regime  are  in  excellent  agree- 
ment with  those  found  from  the  linear  dispersion  relation.  The  exponen- 
tially growing  wave  attains  a sufficiently  high  amplitude  by  X = 350  to 
start  nonlinearly  affecting  the  particle  dynamics  and  the  growth  rate 
begins  to  decrease  until  it  reaches  zero  at  T ~ 420  when  saturation 
occurs . 

The  nonlinear  change  in  the  growth  rate  is  accompanied  by  a less 
pronounced  change  in  the  frequency  shift  prior  to  saturation.  We  find 
that  for  the  parameters  used  in  this  example  at  saturation 

+ a(x)  - ~ 0 in  good  agreement  with  the  qualitative  argu- 

ment in  Section  VI.  The  normalized  field  amplitude, 

E#(x)  * 'e|  EQ(x)/(moc  Wq) , and  the  average  perpendicular  particle  energy 
<Y^>  are  shown  as  functions  of  x in  Fig.  15. 

We  have  found  that  the  energy  efficiency  curves  for  a given  set  of 
parameters  always  display  a maximum.  For  example,  with  u = 0.1,  the 
efficiency  curve  shown  in  Fig.  16  has  a peak  of  40  percent  when 
YQi  = 1.1.  The  peak  results  from  the  competition  between  the  two  satu- 
ration mechanisms  described  in  Section  VI.  As  can  be  seen  in  Fig.  16, 

the  efficiency  curve  has  a cutoff  at  y = Y . = 1.0185.  The  curves 

J.  A.crit 

of  the  linear  growth  rates  and  frequency  shifts,  obtained  by  numerically 
solving  the  linear  dispersion  relation  of  Eq.  (1)  are  shown  in  Fig.  17 
as  a function  of  Y . The  linear  growth  rate  curve  shows  a threshold  at 

Oi. 

y m 1.0185,  identical  to  that  in  Fig.  17.  Figure  17  also  shows  a 
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comparison  between  the  results  of  linear  theory,  Eq.  (1),  and  of  the 

particle  simulations  in  the  linear  regime.  Clearly,  the  linear  growth 

rates  and  frequency  shifts  are  in  excellent  agreement.  For  beams  that 

have  initial  energies  characterized  by  y we  find  that  the 

or  j.,crit 

excited  wave  extracts  free  energy  from  the  beam  until  <Y  > = Y 

a j.,crit 

and  the  system  becomes  linearly  stable  without  particle  trapping  taking 
place.  The  simulations  confirm  this,  since  the  velocity  phase  space 
plots  show  no  evidence  of  trapping  at  the  time  of  saturation  when 
Y0i  ^Yx  Crit’  c^*-s  re8:i-me  a rough  expression  for  the  average  change 

in  y at  saturation  is 

' i 

< &Y  > = v - v 

S ' ox.  Tx,crit 

and  implies  an  efficiency  monotomically  increasing  withyoi;  this  is 

shown  as  a dashed  curve  on  Fig.  16.  In  the  other  regime  where 

Y >>  y • the  saturation  mechanism  is  dominated  by  particle  trap- 
o x i.,crit 

ping  and  we  find  that  <Aya  > g = 2f  yqx  Aw/u^.  This  implies  an  efficiency 

curve  monotomically  decreasing  with  y as  shown  in  Fig.  16  with  f = 1. 

The  two  saturation  mechanisms  described  here  predict  an  efficiency 

maximum;  this  is  also  verified  by  the  results  of  the  particle  simulations. 

Our  saturation  arguments  show  qualitative  agreement  with  the  actual 

efficiency  obtained  numerically.  In  the  YQ±>>  Y ± crit  re8ime  we  have 

arbitrarily  set  f = 1,  implying  that  all  the  particles  are  deeply 

trapped  at  saturation.  From  Fig.  16  we  see  that  in  the  intermediate 

regime  f is  closer  to  1/2.  Similar  arguments  hold  in  the  y ^ y 

° oi  a, crit 

regime.  The  dashed  curves  in  Fig.  16  show  the  competition  between  two 


mechanisms  leading  to  a maximum  in  the  efficiency  curves. 
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Ke  now  present  the  results  of  our  nonlinear  analysis  for  a wide 
range  of  parameters.  In  Fig.  18,  energy  efficiency  as  a function  of 
y is  shown  for  various  values  of  u , u = 0.025,  0.05,  0.15  and  0.35. 

Oi 

In  this  figure  the  initial  wave  frequency  ui^  was  set  equal  to  ck^  with 

n m 1.  Furthermore,  the  initial  relativistic  cyclotron  frequency, 

ft  /y  was  set  equal  to  ck  . These  initial  conditions  on  u , ck  and 
o ox  n on 

ft  /y  were  selected  because  they  occur  near  the  maximum  linear  growth 

o Ox 

rate.  All  of  the  curves  in  Fig.  18  display  a maximum  in  efficiency  in 
the  beam  frame  as  a function  y . This  characteristic  of  the  curves  has 

Ox 

already  been  explained.  A further  feature  of  the  curves  is  that  as  y 


increases,  y . increases  in  agreement  with  the  definition  of  8 

x.crit  ° x.crit 

given  in  Section  III.  From  Fig.  18  we  also  note  that  the  efficiency  in- 


creases at  low  values  of  yq  as  y is  decreased.  This  property  of  the 

curves  can  be  explained  by  noting  that  as  y decreases,  y decreases, 

x,cnt 

resulting  in  a sharp  increase  in  efficiency  as  a function  of  y near 

ox 

y ± c^iC*  The  maximum  efficiency  seems  to  level  off  at  — 40  percent  for 

low  values  of  y , when  w = ck  = ft  /y 

ox  0 n o ox 

Dramatic  increases  in  the  efficiencv  can  be  realized  in  a number  of 


ways.  For  example,  by  slightly  mistuning  the  relativistic  cyclotron 

frequency  away  from  = ck^,  electron  phase  trapping  can  be  postponed 

and  higher  efficiencies  realized.  In  the  regime  where  electron  trapping 

is  responsible  for  saturation,  we  have  shown  that  the  efficiency  is 

roughly  2f  y (w  + Acj  - ft  /y  )/(Y  - 1).  By  choosing  <*>  > ^ /Y 

we  can  expect  the  efficiency  to  increase.  In  Fig.  19  the  efficiency  is 

shown  as  a function  of  (ft  /y  ) /to  for  various  values  of  y and  > 

o ox  o ox 

The  values  for  y and  yQ^  were  taken  from  the  maximum  efficiency  points 
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in  Fig.  18.  Figure  19  shows  clearly  that  higher  efficiencies  can  be 

realized  for  (SI  /y  )/^  less  than  unity.  Since  sufficiently  small 
o ox  o 

values  of  (ft  /y  )/u  lie  outside  the  growth  rate  spectrum,  the 

O OX  O or. 

efficiency  drops  to  zero.  By  mistuning  the  magnetic  field,  the 

efficiency  was  increased  from  40  percent  to  63  percent  in  the  beam  frame 

for  the  case  p = 0.1  and  y = 1.155.  . 

Ox 

A 

The  efficiency  can  also  be  increased  varying  the  external  magnetic 

field  in  time  in  such  a way  that  w - ftQ  (t)/< y^ (t) > is  held  fixed. 

This  procedure  is  similar  to  the  mistuning  approach  discussed,  in  that 

electron  phase  trapping  is  postponed.  Choosing  u = 0.1,  y =1.1  and 

ox 

changing  the  magnetic  field  in  time  such  that  (u  - ft  (x)/<y  (t)  > ) 

o x 

is  fixed  at  its  initial  value  of  0.034,  an  efficiency  of  =s  70  percent 

was  realized.  This  compares  with  an  efficiency  of  ^ 40  percent  when 

the  magnetic  field  was  held  fixed  at  ft  /y  = w = ck  . 

o ox  o n 
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Section  VIII.  Discussion  and  Example 


We  have  shown  that  the  electron  cyclotron  maser  mechanism  can  be 
an  efficient  method  of  producing  short  wavelength  radiation.  The  mecha- 
nism seems  particularly  attractive  for  the  generation  of  submillimeter 
radiation  when  higher  cyclotron  harmonics  are  considered.  Preliminary 
results  at  higher  cyclotron  harmonics  and  also  higher  waveguide  modes 
indicate  that  the  efficiency  remains  relatively  high.  These  results, 
on  the  higher  harmonics,  will  be  published  in  a forthcoming  paper. 

Experiments  on  the  cyclotron  maser  using  intense  relativistic 

electron  beams’-0-*3  have  produced  efficiencies  of  less  than  2 percent. 

The  low  experimental  efficiencies  are  in  good  agreement  with  our 

theoretical  results  and  are  the  result  of:  1)  the  transformation  of 

efficiency  factor  from  the  beam  frame  to  the  laboratory  frame,  that  is 

- 1)/(yo’  (yq  - 1)),  see  Eq.  (B 7 ) , was  typically  small,  < 0.3; 

2)  the  value  of  y in  the  beam  frame  was  not  optimized  for  maximum 
ox 

efficiency,  see  Fig.  18;  3)  the  method  used  to  produce  transverse  energy 
also  introduced  temperature  in  the  beam  distribution,  see  Eqs.  (28), 

(31)  and  (32).  The  combination  of  these  effects  led  to  the  low  observed 
efficiencies.  With  our  improved  understanding  of  the  saturation  levels 
and  mechanisms,  efficiencies  can  be  substantially  improved. 

As  an  illustration  of  the  beam  and  waveguide  parameters  needed  to 
generate  electromagnetic  radiation  at  1.0  cm,  we  consider  the  following 
detailed  example.  Referring  to  Fig.  2,  we  choose  a = 0.25  cm  and 
L “ 2*(2a)  « 3.14  cm.  The  waveguide  width,  L,  corresponds  to  the  cir- 
cumference of  a cylindrical  guide  having  a radius  equal  to  2a.  In  the 
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laboratory  frame  we  have  a cold  electron  beam  with  a total  energy  of 

V = 300  keV  per  particle.  The  electron  velocity  is  partitioned  such 
that  v = 0.69  c and  v - 0.36  c.  Therefore,  the  total  gamma  of  the 

Oil  Ol 

beam  electrons  is  y * (1  - (v  /c)‘  - (v  /c)2)-1-/2  = 1.59  and 
o on  oi 

v = (i  _ (v  /c)2)-1^‘  ■ 1.37.  With  a magnetic  field  B = 11.4  kG, 
on  on  o 

the  electron  Larmor  radius  is  r,  s v / ( I e | B / y m c)  = 0.340a  = 0.085 

L oi  o o o 

cm.  Taking  the  axial  current  to  be  I ■ 860A,  the  beam  thickness  is 

equal  to  twice  the  Larmor  radius  and  beam  width  equal  to  L,  we  obtain 

a beam  density  in  the  laboratory  frame  of  n^  = I/(|e|  vq||(2  r^L))  = 

2.6  * 1011  cm-3.  The  electron  surface  density  is  given  by 

o - n (2  r ) = A. 42  x 1010  cm-2, 
o a L 

Now  in  the  beam  frame,  v'  * 0,  hence,  y'  = (1  - (v1  /c)2)-1/‘  » 

on  Ox  Ox 

Y /y  “ 1.155  and  v'  = (1  - (y*  )-2)1/,‘-  = 0.5  c.  The  electron  surface 

o o II  ox  ox 

density  in  the  beam  frame  becomes  o'  - o y'  /y  = 3.2  x 1010  cm-2. 

o o ox  o ’ 

which  determines  the  modified  plasma  frequency,  u'  = (4  tt  I e 1 2 o'/( m a))1 

b o o 

“ 2.02  x 1010  sec-1.  To  obtain  the  efficiency  in  the  beam  frame,  q', 

we  need  i'  as  well  as  (0  /y'  )/(ck  ).  To  continue,  we  find  that 

o ox  n 

* lel  B0/ (Yq.  mQc)  » 1.75  x 1011  sec-1,  ckn  * 1.88  x 1011  sec-1 

and  (0  /y'  )/(ck)  *=  0.93.  Now  p ' * (ui'/t^y'  ) /cu ' = 0.1,  since  ui’  is  set 
o ox  n b Ox  o o 

equal  to  the  cut-off  frequency,  ck^.  From  Fig.  14  we  find  that  for 

u'  • 0.1,  y'  m 1.155  and  (0  /y'  )/w'  = 0.93  the  efficiency  in  the  beam 
ox  o ox  o J 

frame  is  q'  « 63  percent.  Using  Eq.  (B7)  of  Appendix  B,  the  efficiency 
in  the  laboratory  frame  is  n » n'  y <y ' - l)/(y'  (y  - 1))  = 22  percent 

O Ox  o o 

Since  the  beam  power  in  the  laboratory  frame  is  P - IV  = 258  MW,  the 

b 

radiated  power  becomes  P * n P,  “ 57  MW.  The  spread  in  y'  in  the  beam 

r b x 

frame  due  to  the  self  electrostatic  fields  can  be  estimated  from 


Eq.  (32)  of  the  text.  We  find  from  Eq.  (32)  that  6y'  ,,/y'  = 0.0036. 

j.,self  oi 

From  Eq.  (25)  we  see  that  the  initial  spread  in  Y ' can  be  neglected  if 
it  satisfies  the  following  requirement,  6y'/y’  <<  2 n'(Y'  - 1)y'  = 

± OJ.  OJ.  OJ. 

0.17.  Clearly,  in  this  example,  the  spread  in  y^  due  to  self  fields 
can  be  neglected. 

The  possibility  of  high  efficiency  and  power  levels  at  millimeter 
and  submillimeter  wavelengths  with  the  electron  cyclotron  maser  make 
it  an  attractive  device  for  plasma  heating.  The  frequencies  required 
for  heating  at  either  the  electron  cyclotron  or  upper  hybrid  frequency 
are  in  excess  of  120  Ghz  (X  = 2.5  mm)  for  reactor  type  Tokomaks.23 
The  cyclotron  maser  appears  to  be  the  only  viable  source  for  producing 
the  required  power  efficiently  at  such  high  frequencies.  Plasma  heating 
using  cyclotron  maser  devices  has  been  demonstrated  in  the  USSR  on  the 
TM-3  Toknmak.23  In  these  experiments  the  electron  cyclotron  resonance 
at  X * 5 mm  was  used.  It  has  teen  estimated  that  an  experimental  power 
reactor  would  require  ~~ 100  to  250  MJ  of  supplemental  heating  in  2 to 
6 seconds.30  This  is  a continuous  wave  requirement  of  ^ 50  MW  at  a 
frequency  of  280  Ghz  (X  = 1.07  mm),  assuming  electron  cyclotron  heating 
and  the  reactor  parameters  of  Reference  30.  On  a large  device  these 
requirements  may  be  achieved  by  employing  several  maser  devices  simul- 


taneously. 
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APPENDIX  A 


The  linear  dispersion  relation  can  be  recovered  from  the  self- 

consistent  nonlinear  equations  derived  in  Section  III,  Eqs.  (1)  and  (13). 

To  facilitate  the  calculation,  we  take  the  limit  where  the  Larmor  radius 

of  the  electrons  is  small  compared  to  the  characteristic  transverse  scale 

length  over  which  the  fields  vary,  nx  <<  a.  In  this  approximation  the 

fields  take  the  simplified  form  E = E (t)  cos  (u)  t + a(t))  and  B = 0, 

y o o z 

where  the  x dependence  has  been  removed  because  in  the  small  Larmor 
radius  limit,  k (a  - X)  ~ nx/2  . In  the  linear  approximation , the  magni- 
tude  of  the  particle  velocity,  B^(t),  and  rotational  frequency,  <J/(t), 
are  expanded  as: 


6 (t)  = 6 + B(l)(t), 

X OX  1 

ij>(t)  * n h + (ai) 

o OX 

where  <<  g and  yC1)  <<  Q /v  are  slowly  varying  functions  of 

x ox  o ox 

time  and  are  linear  in  the  wave  amplitude,  E (t) . In  this  regime  it  is 

o 

appropriate  to  set  E (t)  * e exp  (T, t)  and  a(t)  = Aw.  t where  T.  <<  w 

o o L L L o 

is  the  linear  growth  rate  and  Aw  <<  w is  the  linear  frequency  shift. 

Li  O 

Substituting  (Al)  into  the  orbit  equations  (13),  we  obtain  the  following 
slow  time  scale  equations: 

t 


B(l)(t) 

X 


' e ‘ e0  e 


L 


- . — sin  (6wt  + $ ) , 

2 m c o 


ox  o 


,0)(t)  - - 8 y2  B(l>(t) 


OX 


ox  ox  X 


■ , rLC 

1 e I e p 

1 o 

2 6 Y m c 

OX  ox  O 


(A2) 


cos  (5wt  + $ ) (A3) 

o 


where 


w - /y  + Awt  . In  Eq.  (A3)  the  approximation 
o o ox  L 
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Y (t)  ~ v +6  v 3 80)(t)  is  used.  The  linearized  expression 

A OA  OA  OA  A 

for  the  y component  of  particle  velocity  is: 


v^(t)  = c 3 


+ c 


5 sin  / — — t + $ ^ + c B^(t)  sin  ( — — t + q ) 

oi  Vox  7 " \r0A  7 

6 f ^ 1 ^ ( t ' ) dt'  cos  ( — — t + $ ) . 

OA  J V'ni  °/ 


Solving  Equations  (A2)  and  (A3)  for  times  much  longer  than  a growth 


time,  t >>  I/7  yields 


v (t)  = c 8 sin 
y oa 


infe  ‘ + *•) 


i«i  v L / 

4 y m D*-  ) \ 

oj.  o ( ' 

/ (i  + 7 ) 

«<*>  - l 

' ’ n 1 


(1  + y2  ) 


oa  Y 


OJ- 

rL 

- 2 32 

OJ. 

0 

Y^ 

04. 

yoa 

no  (6 

to2  - 

sin 

Yoa  " 

D2 

/ 

: and 

t 

f * 

(1>(t’) 

dt1 

n 6ooF  \ 

o Ll  — 
— 7— Icos  cot 

Y / 

n 1 ’ 


where  to  = 00  + Aco,  , D = 6<o2  + F.  2 

o L L 

o 

equations  governing  the  frequency  shift  and  growth  rate,  Eqs.  (10a) 
(10b),  in  the  linear  approximation  become: 


(w2  - c2  k2  + 2 w Ato,  ) e exp  (T  t) 
o n o L o L 


2 e or. 


10  tu 

o 


t+2it/co  2tt 

/ -7 


d$  v ($  ,t’ ) sin  <0  t' 
o y o 


(A4) 


(A5a) 
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1 e ! a 


\ e0  exp  (:l°  = 


00  W 


t + 2TT/ 


/ -/ 


dl  v ($  ,t')  cos  u t' 
o y O 


(A5b) 


where  v ($  ,t)  is  given  by  (A4) . Since  AwT  /u  , r /w  « 1 and  u ^ ck 
y o LoLo  on 

we  approximate  w2  + c2k 2 by  2 go2  on  the  righthand  side  of  Eq.  (A5).  By 
multiplying  (A5b)  by  >/-T,  adding  it  to  (A5a)  and  performing  the  integrals 
over  t'  and  $ we  arrive  at  the  linear  dispersion  relation: 

io2go  / « 2 


(< 


W.GO  / , 6 \ 

j2  _ c2k2)  = _b / 1 _ _o±  Gj | 

n 2 y 1 (u)  — fi  /y  ) 2 (u  - 0 /y  W’ 

OA  \ ooj.  o o a / 


(A6) 


Jhere  u) 


Wo  + + irL‘  This  result  is  identical  to  the  dispersion 


relation  given  by  Eq.  (1)  for  waves  of  zero  wavenumber,  when 


v <<  (aft  /y  )/n. 
oa  o Oi 
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appendix  b 


One  of  the  objectives  of  this  paper  is  to  determine  the  efficiency 
of  electromagnetic  energy  flux  radiated  by  a beam  undergoing  the  electron 
cyclotron  maser  instability.  The  general  law  governing  the  laboratory 
frame  efficiency  of  conversion  of  beam  to  electromagnetic  power  in  terms 
of  beam  frame  energy  efficiency  is  derived. 

The  analysis  in  the  body  of  the  text  has  been  performed  in  a rather 
special  frame  in  which  beam  axial  velocitv,  v , and  wave  group  velocity. 

Oil 

Vg,  both  vanish.  The  justification  for  choosing  this  frame  is  that  the 
growth  rate  maximizes  when  vq||  = v^;  tl.""efore,  by  transforming  to  the 
beam  frame,  the  axial  beam  velocity  and  the  wave  group  velocity  vanish 
simultaneously.  In  the  beam  frame  the  wave  is  at  cutoff  having  an 
infinite  axial  wavelength,  kQ  - 0.  In  what  follows  the  two  frames  are 
distinguished  by  denoting  beam  frame  variables  by  primed  quantities  and 
laboratory  variables  by  unprimed  quantities.  TU«  notation  in  this 
Appendix  should  not  be  confused  with  that  in  the  t mt  where  unprimed 
quantities  refer  to  the  beam  frame. 

In  the  laboratory  frame  the  axial  electromagnetic  energy  flux 
density,  Sz  = (c/4tt)  (E  * 3)  • e^,  is  found  by  a Lorentz  transformation 
of  the  beam  frame  fields  (derivable  from  Eq.  (6))  and  takes  the  form: 

S * -r-  y2  B (E'(t)  cos  (u't  + a')  sin  k (x  - a))  , (Bl) 

z 4 it  on  on  \ o o n / 


where  6 = v /c,  v is  the  axial  beam  velocity  in  the  laboratorv 

o II  o II  O II 

frame  and  y = (1  - 32  Using  Eq.  (Bl),  the  temporal  and 

o II  o II 

spatial  average  of  the  electromagnetic  power  in  the  laboratory  frame. 
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flowing  normal  to  the  cross  sectional  area,  2aL,  is 


P = — y 
f 8^  on 


on 


(E*  )■ 
o 


aL 


(B2) 


where  L is  the  length  of  the  waveguide  as  shown  in  Fig.  2.  From  con- 
servation of  field  and  particle  energy,  Eq.  (16),  we  obtain  an  expres- 
sion for  the  maximum  electric  field,  occurring  at  saturation,  in  the 
beam  frame: 


. r-  / \l/2  / \ 1 / 2 / 

E ’ = H B U'  L'/Cil'/Y'  ) (Y  - <Y ' > y' 

o,s  O \ O O OJ . ) \ OX  - I \ ox 


-i  r 


(B3) 


Expressing  the  electric  field  amplitude  at  saturation,  Eq.  (33) , in 

terms  of  the  beam  frame  energy  efficiencv,  n'  = (y'  - yy ' l/(y'  - 1) 

\ cx  x s ' ox 

and  substituting  it  into  Eq.  (B2)  gives  for  the  radiated  power  in  the 
laboratory  frame: 


4ir  ‘on 


B“  (p’) 


m2 


n'/Y' 

o ox 


(Y'  - 1) 

ox 


ox 


aL  o', 


(B4) 


The  efficiency  of  conversion  of  beam  to  electromagnetic  power  in  the 
laboratory  frame  is: 


n = v /n 

n V*b> 


where  “ cq  L c 3q||  (yq  - 1)  mQc^  is  the  beam  power,  3 is  the  particle 

-1/2 

surface  density  and  Yq  * (1  - aQ(|  - 3^)  * Combining  Eqs.  (B4)  and 

(B5)  and  using  the  relationship  a(B  p')  2 u)'/(4nm  c2f!')  = a ' 

o o o o o’ 

we  obtain 


o’  (Y'  - 1) 

n » -°  y2  _2i 

a on  (y  - 1) 
o o 


(B6) 
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Applving  the  relativistic  transformations  o = o'  y /y'.  and  y = y /y' 

° o o o or  on  o 01. 

to  Eq.  (B6),  the  power  efficiency  takes  the  simple  form: 


’o  (C  - » , 

° " Yl  <Y.  " 1)  n ' 

0i  o 


(B7) 


For  a more  direct  derivation  of  Eq.  (B7)  we  consider  the  relativis- 
tic Lorentz  transformation  of  total  field  energy  in  the  beam  frame  to 
the  laboratory  frame.  Since  the  field  momentum  in  the  beam  frame  is 
zero,  we  have: 


W,  = y 
f on 


w;  = y n' 
f on 


W 

V 


(B8) 


where  W,  amd  W,  are  the  total  field  and  particle  energy  in  the  labora- 
t b 

tory  frame.  The  power  in  the  electromagnetic  field  can  be  written  as: 


?c  = v W /Az  = v r,'  W'/Az,  (B9) 

f on  f on  b 

where  Az  is  the  axial  extent  of  the  volume  element  which  contains  the 
total  field  energy  W^.  The  total  beam  energy,  W^,  in  the  beam  frame 
can  be  written  as: 


W’  = nK  mr,c2  A Az'» 

b b 01  o 


(BIO) 


where  A = 2 r^  L is  the  cross  sectional  area  of  the  beam  and  Az'  is  the 

axial  extent  of  the  volume  element  in  the  beam  frame.  Substituting 

Eq.  (BIO)  into  (B9)  and  using  the  relationships  Az'  = q,  Az  and 

b b 

Y = y /y ' yields  the  result: 
on  o Oi  y 


Y (Y ' - 1) 

o oi 


f Y'  (y  - 1)  ' b’ 
oj.  o 

which  defines  the  efficiency  transformation  in  Eq.  (B7). 
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Fig.  1 — The  mechanism  responsible  for  the  electron  cyclotron  maser 
illustrated  by  orbits  of  test  particles  in  velocity  space  in  the  presence 
of  a small  external  field,  (a)  initial  particle  positions,  and  (b)  bunched 
particles  after  several  cycles. 
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Fig.  2 — The  electron  cyclotron  maser  configuration 
in  plane  geometry 
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ck0/ckn 

Fig.  4 — Typical  dispersion  curve  for  cyclotron  maser  instability.  This  case 
is  for  T0i  = 1-2,  Wb/v/701  = °-05  wo*  wo  = ckn  = cno/70l  and  C = n = 1 
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(Q0/y01)/(ckni 


(ck0)/(cknl 


Fig.  5 — Contours  of  the  growth  rate  for  the  cyclotron  maser  instability  vs 
the  parallel  wave  number  and  the  cyclotron  frequency.  This  case  is  for  7oi 
= 1.2  and  - 0.05  u>0,  where  coQ  = ckn  and  t = n = 1. 
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Fig.  6 — Particle  trajectories  in  velocity  phase  space  as  found  from 
the  constant  of  the  motion.  The  figure  shows  the  boundaries  of 
the  regions  accessible  to  particles  initially  uniformly  distributed  be- 
tween 0 < X < 27t  with  u i = u . . 


Fig.  7 — The  deformation  of  a monoenergetic  beam  in  the  presence  of 
a constant  amplitude  field  showing  the  actual  particle  positions  in  phase 
space.  The  curve  for  t = tQ  indicates  the  initial  positions  and  the  one 
for  t = tj  the  positions  when  the  beam  indicates  trapping. 
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Fig.  8 — The  particle  positions  in  phase  space  for  a monoenergetic  beam 
deformed  by  a constant  amplitude  electric  field  at  t = t2  when  the  par- 
ticles have  lost  the  maximum  energy  and  are  in  a state  corresponding  to 
saturation. 
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Fig.  10  — The  particle  positions  in  phase  space  from  a simulation  of 
the  cyclotron  maser  instability  in  the  case  of  saturation  by  energy 
depletion,  (a)  t = tj  linear  phase,  (b)  t = t2  nonlinear  phase,  and 
(c)  t = t3  at  saturation. 
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. 11  — The  distribution  function  of  the  electrons  from  a case 
saturation  by  energy  depletion  at  various  times  as  in  Fig.  10 


Fig.  14  — The  nonlinear  phase  shift  and  growth  rate  as  a 
function  of  time  for  a typical  simulation 
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ie  cyclotron  maser  as  a function  of  eneroy 
lisms  of  saturation  and  from  simulations 
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Fig.  17  — A comparison  of  the  growth  rate  and  phase  shift  predicted 
by  linear  theory  and  the  results  from  the  linear  portions  of  the  com- 
puter simulations 
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Fig.  18  — The  beam  frame  efficiency  of  the  electron  cyclotron  maser 
as  a function  of  energy  for  various  densities  when  ckn  = no/70| 
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Fig.  19  - The  beam  frame  efficiency  of  the  electron  cyclotron 
maser  when  the  cyclotron  frequency  is  varied 
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